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Given a real-valued phase-space function, it is a nontrivial task to determine whether it corre- 
sponds to a Wigner distribution for a physically acceptable quantum state. This topic has been 
of fundamental interest for long, and in a modern application, it can be related to the problem 
of entanglement detection for multi-mode cases. In this paper, we present a hierarchy of com- 
plete conditions for a physically realizable Wigner distribution. Our derivation is based on the 
normally-ordered expansion, in terms of annihilation and creation operators, of the quasi-density 
operator corresponding to the phase-space function in question. As a by-product, it is shown that 
the phase-space distributions with elliptical symmetry can be readily diagonalized in our represen- 
tation, facilitating the test of physical realizability. We also illustrate how the current formulation 
can be connected to the detection of bipartite entanglement for continuous variables. 

PACS numbers: 03.65.Ta, 42.50.-p, 03.65.Ud 



Oh' 



> 
(N 

O 

OO 

O 



X 



I. INTRODUCTION 

The phase-space formulation of quantum mechanics 
has continued to be of great interest ever since its first 
introduction by E. Wigner Remarkably, J. E. Moyal 
further elaborated on the so-called Wigner- Weyl corre- 
spondence through an attempt to interpret quantum me- 
chanics as a general statistical dynamics 0- The phase- 
space treatment of quantum mechanics is particularly in- 
teresting in that it provides a valuable insight into the 
issue of quantum-classical correspondence Q . It exhibits 
both the similarities and the differences between the 
quantum and the classical descriptions of physical world. 
Although the original scheme was intended to shed light 
on quantum mechanics from a classical perspective, the 
inverse problem, i.e., the description of classical mechan- 
ics in Hilbert space from a quanturnperspective, has also 
attracted much interest [1, |3, H, S Q ■ 

The most prominent difference between the quantum 
and the classical theories may be characterized by the 
uncertainty principle. According to it, every classically 
admissible distribution in phase space is not always al- 
lowed by quantum mechanics. In classical mechanics, 
any positive-definite real- valued functions, which may be 
interpreted as joint-probability densities for having defi- 
nite values of both the position and the momentum, are 
possible in phase space. On the other hand, in quantum 
mechanics, due to the canonical commutation relation 
[x,p] = ih, the uncertainty relation AxAp > ^ must be 
at least satisfied by the distributions. 

However, it is also known that the single uncertainty 
condition AxAp > | alone does not guarantee that the 
phase-space function corresponds to a quantum mechani- 
cally acceptable state [l^ llH . In order to more deeply 
understand how the quantum principles are reflected in 
phase space, it may be of crucial importance to have com- 
plete conditions, desirably in analytic form, by which to 
determine whether a given real-valued phase-space func- 
tion may describe a legitimate quantum state or not. 
This actually defines a substantially nontrivial task in 



contrast to the inverse problem: When a certain quan- 
tum state described by a density operator p is given, it is 
now a well-known, standardized, procedure to obtain the 
phase-space distributions corresponding to various oper- 
ator orderings jU]. However, even the set of classical 
phase-space distributions that are allowe d by quantum 
mechanics is not yet completely identified [13j . 

In this paper, we derive a hierarchy of complete con- 
ditions for a legitimate Wigner distribution of d degrees 
of freedom. It is done by using the correspondence be- 
tween a given real-valued function in phase space and a 
Hermitian operator pq in Hilbert space, which may be 
called a quasi-density operator [(||. In particular, the 
quasi-density operator pq is constructed as a power-series 
expansion of annihilation/creation operators in normal- 
ordering. We then require that the quantum fidelity of 
the quasi-density operator with an arbitrary pure state 
must be non-negative as a sufficient and necessary con- 
dition for the positivity of pq. We compare these condi- 
tions with the Kastler-Loupias-Miracle-Sole (KLM) con- 
ditions [13, [3, the only systematic conditions previously 
known to our best knowledge, and remark that some 
practical advantages may arise from our novel approach. 
As a by-product, we show that the quasi-density operator 
corresponding to a phase-space distribution with ellipti- 
cal symmetry can be always diagonalized in the general- 
ized Fock-state basis, which may facilitate the test of pos- 
itivity. We also illustrate the utility of the derived con- 
ditions by examining the Wigner function corresponding 
to a two-mode entangled state under partial transpose 
flG^I , which may be interpreted as an observable criterion 
for bipartite entanglement. 

This paper is organized as follows. In Sec. II, some 
preliminary facts are briefly introduced to be used later, 
and in Sec. HI, the previously known KLM conditions 
are reviewed in a heuristic manner. In Sec. IV, a set 
of complete conditions for legitimate Wigner distribu- 
tions is derived based on the normally-ordered expansion 
of the quasi-density operator, and the case of elliptical 
symmetry in phase space is more speciflcally considered. 
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These results are further extended particularly to two- 
mode cases in Sec. V, and the formalism is illustrated by 
an example relating to the problem of detecting bipartite 
entanglement. In Sec. VI, the results are summarized 
with concluding remarks. 

Before going to the main part, we need to clarify our 
usage of operators in this paper. The annihilation and 
the creation operators, a and a\ for one degree of free- 
dom satisfy the well-known boson commutation relation, 
[a, a^] = 1. When the annihilation operator is decom- 
posed into real and imaginary parts as a = X + lY , the 
two Hermitian operators X = and Y — , which 
are known as the quadrature amplitudes in quantum op- 
tics, satisfy i^] = |. In other words, setting fi. = 1, 
the quadrature amplitudes can be related to the position 
and the momentum operators as i = V2Xandp = V2f. 
In this paper, we will take into consideration only the 
quadrature operators X and Y and their corresponding 
values in phase space, not the canonical operators x and 
p. However, the translation of the results into terms re- 
lated to X and p is rather straightforward by considering 
the numerical factor \/2. 

II. PRELIMINARIES 

Let us first consider the case with one degree of free- 
dom. It is well known that an arbitrary bounded op- 
erator F with a finite Hilbert-Schmidt norm, = 
Tr{_F^^i^} < CX3, can be represented in an integral form as 

F^^ j dP\b\\)Cp{\l (1) 

where b{\) = e'^"^"'^ ° is the displacement operator. 
The complex- valued function Cf(A) = Tt:{FD{\)} is 
usually termed the characteristic function of the oper- 
ator F If we further proceed with the normal- 
ordered expansion of the displacement operator using the 
Baker-Campbell-Hausdorff relation, b{X) = e^'^^"^'° = 

e e e 2 ^ the operator F can be cast into the form 
F^iVC^^at^a", (2) 

TT 

where the coefficients C„i„ are given by 

C™„ = / rf2AA'"A*"e' — Cf(A). (3) 

mW. J 

When the operator F is Hermitian, F — F\ then the 
characteristic function satisfies the property Cp{\) = 
Cp{—\), according to which Cnm = Cmn follows. It 
is also known that the coefficients Cmn are all finite, 

bounded as | Cmn \ < ' that the normally- 

ordered power series in Eq. ^ converges well to the op- 
erator F [l3|. 



Specifically, when the bounded operator F is a quan- 
tum density operator p, the characteristic function 
C(A) = Tr{/5Z)(A)} is used to define the Wigner distribu- 
tion T/F(q:) via a Fourier transform in phase space, i.e., 

Wia) = J d2Ae"^'^"*-^C(A). (4) 

The quantum state of any one-dimensional system can 
be thus represented either by the operator p or by the 
c-number Wigner function W{a), as the correspondence 
in Eq. (jj]) is one-to-one. For a legitimate quantum state, 
the density operator p is positive semidefinite, i.e., all its 
eigenvalues are nonnegative, along with the trace condi- 
tion Tr{/5} — 1, which is our concern in this paper. 

legitimate Wigner distribution— Suppose now that 
a certain real-valued function Wt{a) — Wt{ax,ay) is 
given to be tested whether it corresponds to a legitimate 
quantum state, or more precisely, a legitimate Wigner 
distribution. One can first take its complex Fourier trans- 
form to obtain the characteristic function as 

Ct{X) - J d^ae^^'-^'^Wtia). (5) 

Now, from Eq. ([1]), the quasi-density operator pq corre- 
sponding to the given Wt{a) is constructed as 

P^ = l^J d^^b\X)Ct{X). (6) 

The question at hand is to check if pg is positive semidef- 
inite with the trace condition Tr{pq} = 1. The trace 
condition is rather simple to test, as Tr{pq} = Ct(0) = 
/ (PaWt{a) = 1, i.e., it requires that Wt{a) is integrated 
to unity over the entire phase space. Thus, we will focus 
only on the positive semidefiniteness of fSq throughout 
this paper. 

Note that other types of phase space distributions can 
be treated in a similar way. Supposed one is given a real- 
valued function Wf{a) = Wf{ax, cty) and asked whether 
it corresponds to a specific s-parametrizcd distribution in 
phase space 12] . The Wigner distribution corresponds to 
the case of s = 0, and other important distributions in 
quantum optics are the Glauber-P function for s = 1 and 
the Q-function for s = — 1 . One would then calculate the 
characteristic function as 

Ct(A) = e-5l^l' J d2ae^"*-^*"W^/(a), (7) 

to obtain the quasi-density operator in the form of 
Eq. ([5]). In this paper, we will deal only with the case of 
s — 0, that is, the Wigner distribution. 

III. KLM CONDITIONS 

To begin with, let us briefly address the KLM con- 
ditions for the positivity of the quasi-density operator 
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Pq [13 ■ If Pq is positive semidefinite, the condition 
Tr{/ vPg} ^ must be fulfilled for an arbitrary operator 
/. Let us particularly take / as a discrete sum of the dis- 
placement operators, i.e., / = Yll^^i where Ai 
and are arbitrary complex numbers. Then, using the 
relation 



D{a)D{l3) = D{a + I3)e 



i(a/3*-a-/3) 



(8) 



and Tr{D{a)} = 7r(5^(a), along with Eq. it follows 



n 

Tr{/VVj = - E > 0, 



where 
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(9) 



(10) 



As the inequality Eq. Q must be fulfilled for arbitrary 
Ai^s, every n x n matrix {My} (n = 1,2,...) must be 
positive semidefinite, which becomes the necessary con- 
dition for pg to be positive semidefinite. 

According to KLM [l3, [l^, the converse is also true, 
i.e., the positive semidefiniteness of all matrices {Mij} 
becomes sufficient for that of pq. In our heuristic argu- 
ment, this sufficiency may be seen, with less mathemati- 
cal rigor, as follows. As implied by Eq. ((T]), an arbitrary 
operator / can be expressed as a sum, more precisely 
an integral, of the displacement operators. Then, fol- 
lowing the same procedures as before, the requirement 
of Tr{/^/pg} > for every operator / leads to a simi- 
lar condition to Eq. ([9]), where the discrete sum is only 
replaced by a continuous integration. 

The KLM conditions may be understood as a quantum 
version of the Bochner theorem [3, [3 ■ The " classical" 
Bochner theorem states that every positive-definite func- 
tion is a Fourier transform of a positive finite Borel mea- 
sure. In our problem, it implies that if the phase-space 
function is a classical probability density {Wt{a) > 0), 
then the matrix {M-j = Ct{£,i — ^j)}, instead of the one 
in Eq. (|10p . must be positive semidefinite. Thus, the only 
difference between quantum and classical cases is the ad- 
ditional factor e^^i^^i of Eq. pUj) in quantum cases, 
which obviously arises due to the commutation relation 
[a,a^ = 1 through Eq. 

On constructing the n x n matrix {Mij}, one has to 
show that Mij > for every choice of {i = 1, • • • , n) 
to confirm the positivity of pq. In many cases, the test 
may thus amount to the optimization problem involving 
2(n — 1) real independent variables corresponding to the 
complex variables ^i — ^i+i (i = 1, • • • , n — 1), which be- 
comes increasingly hard for a growing number of n. We 
will now derive another set of complete conditions for the 
positivity of pq that can avoid such a problem. Further- 
more, our approach will turn out to make it possible to 
directly evaluate the eigenvalues of Pq rather easily for 
the cases of elliptical symmetry in phase space. 



IV. COMPLETE CONDITIONS FOR WIGNER 
FUNCTION 

Given the phase-space distribution Wt{a), or equiva- 
lently, the characteristic function Ct(A) in Eq. one 
may represent the quasi-density operator in a normally- 
ordered form given by Eq. ([2]) as 
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771, n 



.Am n 
7nn^ u . 



(11) 



Here, the coefficients are given by 

Cmn = ^V / rf'AA™A*"e-^a(A), (12) 



or, alternatively in terms of the phase-space function 
Wt{a), by 



27r(-2)' 



[a"-"S'(m,m-7i,2|anVKt(a)^l,3) 



where the kernel S is defined by 



k 

hi 
hi 



{k + hy.kl 
iFi[h + l,h + l,-x]. (14) 



(iFi refers to the hypergeometric function.) 

With the identification of the coefficients Cmn for 
a given distribution Wt{a), the next step is to check 
whether the quasi-density operator in Eq. (jlip is posi- 
tive semidefinite. This can be done using the fact that 
a Hermitian operator H is positive semidefinite if and 
only if > for every pure state |^) S H, where 

H is the Hilbert space under consideration. Express- 
ing an arbitrary pure state in the Fock-state basis as 
I'l') ~ T,kDk\k) {a^a\k) = k\k)), this condition reads 



(15) 



k,k' 



where 



afefc 



1 



-k-\-7n 1 



(16) 



(M = max{0, k — fc'}). Therefore, the positive semidefi- 
niteness of the matrix {a^k'} becomes the sufficient and 
necessary condition for the legitimate Wigner distribu- 
tions. 

As a matter of fact, akk' is simply the matrix el- 
ement of the quasi-density operator in the Fock state 
basis, akk' = {k\pq\k'). In this regard, if one di- 
rectly calculate akw from Eq. ^ using the matrix el- 
ement of the displacement operator, {k\]D^ {\)\k') = 



4 



'|!(_A)'=-'='e-l^lV24'"''^(|A|2), where 4^^(x) 
associated Laguerre polynomial 12], it follows 



IS an 



akk' 




fc-fc' 



-a{x)Lt''\\x\') 



In this paper, instead of Eq. pT)) . we are going to 
deal with the expression in Eq. (|16p . which is more di- 
rectly associated with the normal-ordered operator form, 
Eq. pT|) . by way of the coefficients Cmn- Note that the 
condition {akk'} > has a clear physical interpretation 
as the positivity of quantum fidelity of the quasi-density 
operator with arbitrary pure states. This is a full gen- 
eralization of the approach taken by Manko et al. for a 
particular example in 

Clearly, the condition {akk'} > does not involve the 
optimization problem in contrast to the KLM conditions. 
Instead, all relevant information is incorporated by eval- 
uating the coefficients Cmn of the quasi-density operator 
Pq in the form of Eq. or Eq. ((T3]). The positive 
semidefiniteness of the matrix {akk'} can be character- 
ized only in terms of matrix determinants by Sylvester's 
criterion [l^ . That is, all principal minors constructed 
from the matrix {akk'} must be nonnegative for a given 
distribution Wt{a) to represent a physically realizable 
quantum state. 

Case of elliptical symmetry— Let us now consider 
the case in which the distribution function Wt{ax,ay) 
possesses elliptical symmetry in phase space, more pre- 
cisely, the value of Wt depends only on the parameter 



where a'^ = ax cos 4> — 
ay sin 0, and a'y = ax siw (f> + ay cos. (j). {a,b: real con- 
stants) In other words, the distribution is centered at the 
point {I3x,/3y), and the major/minor axes of the ellipse 
are rotated by an angle (f>. The general elliptical distri- 
bution can be transformed to a standard one centered at 
origin with x- and y-axis as major/minor axes by a uni- 
tary operation as UpqW. Here, U is the displacement 

operator followed by the phase-shift, U = e~*''^°*"D(— 
As the unitary transformation does not change the pos- 
itivity of the Hermitian operator, one may focus on the 
positivity of the elliptic distribution only in a form of 



Wt{ax,ay) = Wtiy + without loss of generality. 
Furthermore, let us now consider a rescaled bosonic 

-Y, where X and Y are 



operator as a' = y/jX + 

the original quadrature operators [20|. The new oper- 
ator a' and its conjugate a'^ obviously define a bosonic 
mode as the commutation relation [a' , a'^] — 1 holds. 
The corresponding rescaled distribution then becomes 
Wt{ax,ay) = Wt(^^|a|), i.e., it now possesses the 
circular symmetry in phase space. In this case, the 
quasi-density operator becomes diagonal in the Fock- 
state basis, as the coefficients Cmn = ioi m ^ n in 
Eq. ([I2D or in Eq. (O. ( Note that [a^"'a"',a'fa] = 
for every integer m.) When pq is expressed as pq — 
i Cmma'''"'o™, the diagonal terms become akk = 



{k\pq\k) = ^ X]m=o (k-L)'. ^'^"'-' which of course corre- 
spond to the eigenvalues of pq. The test of positivity of 
Pq thus becomes relatively easy for the elliptical distri- 
-3utions using our method. 

( X7J^s an example, let us consider the distribution intro- 
- iuced by Manko et al. in Q , which is obtained by rescal- 
ing the initial Wigner distribution as ax — *■ XxCtx and 
ay —^ Xyay. Mauko et al. particularly showed that for 
\x = Xy — X, the uncertainty relation AxAp > § is 
fulfilled with the condition A < 1, and that the origi- 
nal Wigner distribution for the Fock state however, 
becomes unphysical for very small A. 

In fact, using our method, one can easily check that the 
rescaled Wigner distributions become unphysical for any 
values of A 7^ 1 . The deformed characteristic function for 
11) reads 



(18) 



where A = y/XxXy. Then, the coefficients in Eq. ([T^ are 
obtained as 



C„ 



(-l)™7r 



2A2 



rn+l 



(A^ - 1 - 2m), (19) 



to!(1 + A2) V1 + A2, 
and the eigenvalues of the quasidensity operator given by 
2A2 



^2 /1_A2^'"' 



akk 



(1 + A2)3 V1 + A2 



[4fcA2 - (1 - a2)2] ,(20) 



where k — 0, 1,2,---. With A < 1, the eigenvalue 
ckoo always becomes negative. On the other hand, with 
A > 1, the eigenvalues akk always become negative for 



k > 



(A"-l)" 
4A2 



Therefore, only the trivial rescaling, A = 1, 



preserves the physical realizability of the Wigner distri- 
butions. 



V. MULTI-MODE CASES 

It is rather straightforward to extend the previous re- 
sults to cases of d degrees of freedom. Specifically, for 
two-mode cases, the quasi-density operator may be rep- 
resented as 

TT ^ — ' 
mn 



where 

Cminim2n2 = 



(-1)"' 



milriil 



d^x,xT'XT' 



The characteristic function Cf (Ai, A2) in Eq. (|^^ can be 
obtained in terms of the two-mode phase-space distribu- 
tion Wt{ai, Q!2) as 



Ct(Ai,A2) 



n\ i2 AiQ*-A*c 
\d aiC ' ' ' 



7=1.2 



Wt(ai,a2),(23) 
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which leads to an alternative expression for Cminim2n2 
as 



r' 



n 

i=l,2 



27r(-2)' 



n [d- 

i=l,2 



2 rrti—ni ry^ 



where the kernel S is defined in Eq. (|14p . Taking similar 
steps, one obtains the condition for the positive semidef- 
initeness of the quasi-density operator as that of the ma- 
trix {akk'}, where 



1 



ki 



k2 



mi— m2— A/2 



(fci - mi)!(fc2 - 1712) 



(25) 



{Mi = max{0, fci — fc^},i — 1,2). Note that we have 
used the collective indices k = {fci, ^2} and k' = {fci, ^2} 
in Eq. In constructing the matrix {akk'}; these 

indices may be arranged in the increasings order of the 
sums fc = fci + fc2 and fc' = k[ + k2, respectively. When 
the sum is fixed, e.g., fc = fci + ^2, one takes the order as 
{fci, fc2} = {0, fc}, {1, fc — 1}, • • • , {fc, 0}, and similarly for 
{k'l , k'2}. 

Detection of bipartite entanglement— As an illus- 
tration, let us consider the problem of detecting bipartite 
entanglement for continuous variables (CVs). The entan- 
glement criteria so far known for CVs are all based on the 
partial transposition (PT) [HI [11 [H [13, [H] . When a 
two-mode state p is separable, it is represented by a form 
p — J^iPiPi'' ® P2*^ where the state pj*^ refers to the 
subsystem j = 1, 2. Under partial transposition (PT) for 
the subsystem 2, a separable state still remaiiis positive, 
therefore it describes a certain physical state [iq . 

In our formalism, for a given two-mode entangled state 
with the Wigner distribution lVt(ai,a2), one can take 
the partially transposed distribution, Wt{ai,a2) f2ll.[25j. 
If the state becomes negative under PT, then the distri- 
bution Wt{ai, must violate at least one of the non- 
negative conditions for the matrix {akk'}, which demon- 
strates the presence of entanglement. Note that the coef- 
ficients C^j^^j^^^^ for the partially transposed distribu- 
tion are given by the original coefficients as C^„_^^^„^ = 
Crmmn2ni2 through Eqs. ([22]) and 

For instance, let us consider the two-mode entangled 
state I*) = a|00) + /3|11). The determinant of the 3 x 
3 matrix {akk'} for the partial-transpose distribution is 
obtained in terms of the original coefficients Cmi7iim2n2 
for the state l^*) as 



1 



.100 > 



Coooo 


Cqoio 


Coioo 


Coooi 


C'oooo + Coon Coioi 


Ciooo 


C'lolo 


C'oooo + C 




\ 




(: 


a/3* 


<o. 




a*/3 J 





(26) 



which is negative for any nonzero values of a and (3. 



The above inequality in Eq. ([26)) implies that by mea- 
suring several coefficients of Cjninin2m2 for a given two- 
mode state, one can verify bipartite entanglement. To ex- 

TOiperiif^^t^]|§^]oTzlt^^i^hi!^)^ values, however, it seems that 
one must construct first the two-mode Wigner distribu- 
tion via the honiodyne^|5|yiography 26|, as seen from 
Eq. This may be regarded as a practical disadvan- 

tage, although one could not rule out the possibility to 
construct the coefficients through an alternative, desir- 
ably more efficient, experimental method. On the other 
hand, there is also a certain advantage in our formula- 
tion regarding entanglement detection. Since we formu- 
lated the complete conditions for Wigner distributions 

+™friafti''th4'2j^SS!tivity of the quantum fidelity of the quasi- 
density operator with respect to pure states in the Fock- 
state basis, one can effectively constrain the dimension 
of the matrix {akk'} to take into account. Specifically, 
when the total excitation of the state is bounded by N, 
one just need to construct at most N' x N' matrix, where 
N' = (^+iK^+^) ^ This was actually done for the above 
example of l^") = a|00) -I- /9|11), for which a further re- 
duction of the dimension of the matrix was possible in 
Eq. (1261). 



VI. SUMMARY AND CONCLUDING 
REMARKS 

In summary, we have derived a set of complete condi- 
tions for a legitimate Wigner distribution of d degrees of 
freedom, based on the normally-ordered expansion of the 
quasi-density operator in terms of annihilation/creation 
operators. It was argued that this set may provide a 
practical advantage over the previously known KLM con- 
ditions. Furthermore, it was shown that the phase-space 
distributions with elliptical symmetry can be rather eas- 
ily diagonalized in our formalism, thereby facilitating 
the test of positivity of the quasi-density operator. The 
derivation was specifically extended to two-mode cases 
and it was illustrated how the conditions can be used in 
detecting bipartite entanglement for CVs, along with the 
discussion on its experimental implementation. 

In this paper, we formulated the complete conditions 
for Wigner distributions using the positivity of quantum 
fidelity of the quasi-density operator with arbitrary pure 
states particularly in the Fock-state basis \n). In one per- 
spective, it would be worthwhile to pursue this quantum- 
fidelity-based approach using generalized basis states 
(Sn \bn){bn \ — I), particularly to identify other types of 
phase-space symmetry for which a diagonalization rather 
easily follows. In another perspective, it was recently 
proved that the positive semidefiniteness of a Hermitian 
operator is equivalent to the satisfaction of all uncer- 
tainty relations [13]. It was further shown that a sep- 
arability inequality can be systematically derived for any 
given negative PT entangled state. In this respect, it 
also seems worthwhile to derive another equivalent set of 
conditions for admissible Wigner distributions purely in 
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terms of uncertainty relations in phase space, which is 
left for future work. 
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